arXiv: 1509.00560vl [astro-ph.EP] 2 Sep 2015 


When Comets Get Old: A Synthesis of Comet and Meteor 
Observations of the Low Activity Comet 209P/LINEAR 

Quan-Zhi Ye 

Department of Physics and Astronomy, The University of Western Ontario, London, 

Ontario N6A 3K7, Canada 

qye22@uwo.ca 

Man-To Hui 

Department of Earth, Planetary and Space Sciences, University of California at Los 
Angeles, 595 Charles Young Drive East, Los Angeles, CA 90095-1567, U.S.A. 

Peter G. Brown, Margaret D. Campbell-Brown, Petr Pokorny, and Paul A. Wiegert 
Department of Physics and Astronomy, The University of Western Ontario, London, 

Ontario N6A 3K7, Canada 

Centre for Planetary Science and Exploration, The University of Western Ontario, London, 

Ontario N6A 5B8, Canada 

and 

Xing Gao 

No. 1 Senior High School of Urumqi, Urumqi, Xinjiang, China 
Received_; accepted_ 




ABSTRACT 


It is speculated that some weakly active comets may be transitional objects 
between active and dormant comets. These objects are at a unique stage of the 
evolution of cometary nuclei, as they are still identifiable as active comets, in 
contrast to inactive comets that are observationally indistinguishable from low 
albedo asteroids. In this paper, we present a synthesis of comet and meteor ob¬ 
servations of Jupiter-family comet 209P/LINEAR, one of the most weakly active 
comets recorded to-date. Images taken by the Xingming 0.35-m telescope and 
the Gemini Flamingo-2 camera are modeled by a Monte Carlo dust model, which 
yields a low dust ejection speed (1/10 of that of moderately active comets), dom¬ 
inance of large dust grains, and a low dust production of 0.4 kg • s“^ at 19 d after 
the 2014 perihelion passage. We also hnd a reddish nucleus of 209P/LINEAR 
that is similar to D-type asteroids and most Trojan asteroids. Meteor observa¬ 
tions with the Canadian Meteor Orbit Radar (CMOR), coupled with meteoroid 
stream modeling, suggest a low dust production of the parent over the past few 
hundred orbits, although there are hints of a some temporary increase in activity 
in the 18th century. Dynamical simulations indicate 209P/LINEAR may have 
resided in a stable near-Earth orbit for ~ 10^ yr, which is signihcantly longer than 
typical JFCs. All these lines of evidence imply that 209P/LINEAR as an aging 
comet quietly exhausting its remaining near surface volatiles. We also compare 
209P/LINEAR to other low activity comets, where evidence for a diversity of the 
origin of low activity is seen. 

Subject headings: comets; individual (209P/LINEAR, 252P/LINEAR, 289P/2003 
WY25 (Blanpain), 300P/2005 JQ5 (Catalina)), meteorites, meteors, meteoroids. 
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1. Introduction 


Dormant comets are comets that have depleted their volatiles in the near surface layers 
but may still possess an ice-rich interior. It is not easy to study these objects directly, 
as their optical properties are indistinguishable from those of some of their asteroidal 
counterparts. Dormant comets among the population of near-Earth objects (NEOs) are 
particularly interesting, as they may have a significant contribution to Earth’s history. It 


(e.g. 


Fernandez et ah 

2002; 

DeMeo & Binzel 

2008) 


The dynamical lifetime of common JFCs is about 10® yr flLevison fc Duncan 


1994 . 


The physical lifetime of 


(e.g. 


Di Sisto et ah 


dlometer-sized JFCs, however, is estimated to be only a few 10^ yr 


20091) • It is therefore evident that a typical JFC, presuming it does 


not fragment or split, would spend most of its time as a dormant comet. The details of 
the active-dormancy transition remain nebulous, but classical understanding of cometary 


evolution argues that the transition might include a period of low or inter mitten 


activity, possibly due to the buildup of dust mantles on the surface (c.f. 


Jewitt 


come t ary 


2904 


Hence, it is natural to speculate that some weakly active comets may be active-dormancy 
transitional objects. From an observer’s perspective, these objects are at a unique stage of 
the evolution of cometary nuclei, as they are still observationally identifiable as physical 
comets, as opposed to completely dormant comets that are indistinguishable from low 
albedo asteroids. 


We dehne a low activity comet as a comet where the absolute total magnitude. Mi, 
is higher (fainter) than the absolute magnitude of a dark asteroid (dehned by V-band 
geometric albedo = 0.1) of equivalent effective body (nucleus) size. The physical 
implication of this definition is that the cometary activity is so low, that the comet would 
be recognized as a dark asteroid (pv < 0.1) if extended cometary features are unresolvable 
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to an observer. Mathematically, the dehnition can be expressed as 


Ml > 16.6 — 5 log 


R 


■N 


1 km 


( 1 ) 


where is the effective nnclens radius. Among the 121 comets with constrained nucleus 
sizeji], we hnd 9 comets meeting our definition of low activity comets (Table [T]) of which 8 
are near-Earth JFCs. 

What are the nature and the origins of these comets? To answer this question, we 
need to look at their physical and dynamical properties. In particular, we note four of 
these comets - namely 209P/LINEAR, 252P/LINEAR, 289P/2003 WY25 (Blanpain) and 
300P/2005 JQ5 (Catalina) - can produce meteor showers currently observable at Earth. 
Meteor showers are caused by cometary dusts ejected in past orbits of the parent, therefore 
meteor observations have the potential of enhancing our understanding of the physical 


history of the parent, as de monstrated in t 


of 55P/Tem pel-Tuttle (e.g. 


comets (e.g. 


Yeomans 


Babadzhanov et al 


2012 


1981 


le investigat ion of the present and past activity 


Brown 


199911 and a couple o f potential dormant 


Kokhirova &: Babadzhanov 


201 ,^ . 


In this paper, we focus on one particular comet in our list, 209P/ LINEAR 


209P/LINEAR is am ong the most weakly active comets ever recorded (e.g. 


Ishiguro et al 


(e.g. 


Schleicher 


Jenniskens 


20141: 


2015h and is associated wit h a new meteor shower, the Camelopardalids 


2014J: 


Madiedo et al. 


2014). What makes 209P/LINEAR ideal in studying 


cometary dormancy transition is (1) the close approach to the Earth of the comet during 
its 2014 perihelion passage, reaching ~ 0.05 AU from the Earth where it had brightened to 
1/ ~ 11 magnitude; and (2) the simultaneous encounter of a series of dust trails produced 


^The nucleus sizes of these 121 comets are extracted from the JPL Small-Body Database 


(http://ssd.jpl.nasa.gov/sbdb_query.cgi) on 2015 June 3. 
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by the comet in its past orbits. These two events provide a rare opportunity to look 
at a potential comet-asteroid transitional object from two complementary approaches. 
Therefore, we observe 209P/LINEAR itself (§2) as well as the associated meteor activity 
(§3) to characterize the current state and recent history of the comet’s activity. The 
observations are coupled with the results from numerical simulations to understand the 
nature and origin of 209P/LINEAR (§4). We also discuss the implication of our results to 
the state of to other low activity comets through the examination of 209P/LINEAR. 


2. The Comet 

2.1. Observation 

Imaging observations were conducted with three facilities at three different epochs. 

The observations and reduction procedures are summarized below and tabulated in Table [2l 

1. Gemini North -|- Gemini Multi-Object Spectrographs (GMOS) camera at 2014 April 
9.25 UT. This is a single frame taken as a snapshot observation. The observation 
was conducted relatively early in the active phase of 209P, making it suitable for 
examining the initial activation of the comet. 

2. The 0.35-m telescope -|- QHY-9 camera at Xingming Observatory on 2014 May 18.75 
UT. Around this date, the viewing geometry was favorable for separating dust of 
different sizes and emission epochs. The observation was conducted without hlters 
and was processed using standard procedures (bias and dark frame subtraction, flat 
frame division). 

3. Gemini South -I- Flamingo-2 (F-2) camera on 2014 May 25.94 UT. Around this date, 
the Earth was close to the comet and was near the orbital plane of the comet. The 
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observation was conducted in the Ks band, with 15 s of exposure of each frame. 
The telescope was nodded in the direction perpendicular to the tail axis, to avoid 
contamination from the tail signal at the sky subtraction stage. As the comet was 
moving at a fast rate of ~ 18"/min (or 25 pix per frame), we opted for the non-guided 
non-sidereal tracking mode to avoid frequent changes of guide stars. Because of 
this, a small fraction (< 5%) of frames suffer from poor tracking and are discarded. 
At the end, a total of 41 frames were useful for later analysis. The data reduction 
is performed with the Image Reduction and Analysis Facility (IRAF) supplied by 
Gemini. 


2.2. Results and Analysis 


2.2.1. Start of Cometary Activity and General Morphology 


Previous researches flHergenrother 


2ni4J: 


Ishiguro et ah 


20151) found that the activity of 


209P/LINEAR started at a small activation distance of rn = 1.4 AU. With the CMOS image, 
we conduct an independent check of the start time of activity of 209P /LINEAR. This is 


done by comparing the surface brightness prohle to a synchrone model (IFinson fc Probstein 


19681 ). We estimate the start of activity occurred no later than late February 2014 or a 


lead time of r ~ 50 d, where 209P/LINEAR was at rn = 1.4 AU (Figure [T]). This is in 
agreement with previous results. 


Composite images taken by Xingming 0.35-m telescope and Gemini F-2 on May 18 
and 25 are shown as Figured In the optical image from Xingming, 209P/LINEAR showed 
a symmetric coma measured 6 — 7" (or about 50% larger than mean Full-Width-Half- 
Maximum or FWHM of background stars) in size and a mostly straight dust tail extended 
beyond the held of view. In the near infrared image from F-2, the nucleus, with the same 
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FWHM compared to background stars, is clearly separated from the coma. The coma is 
signihcantly elongated along the Sun-comet axis, with the sunward side extending ~ 5" or 
~ 230 km towards the solar direction. 


2.2.2. Modeling the Dust 


To understand the dust properties, we model t 


dust model evolved from the one used in 


Ye&Hui 


le ob servations using a Monte Carlo 


20141) . The dynamics of the cometary 


dust are determined by two parameters: the ratio between radiation pressure and solar 


gravity, = 5.7 x 10 where pa the bn 


k density of the c 


ust an d Od the diameter 


model proposed by 


1 SI units ( 

the dust. 

Wvatt & WhionlJ 

1959: 

Burns et ah 

1979 

): and the initial 

The latt 

er is : 

bund following the philosophy of the physical 

Crifo & Rodionov 

(1997 

), is dehned as 


Vei = Vol3]l^ cos z ■ V (2) 

where Vq is the mean ejection speed of a dust particle of /drp = 1, x is the local solar zenith 
angle, and v follows a Gaussian probability density function: 


P{u)=U{l,al) (3) 

where is the standard deviation of v. The T’(zz) function heuristically accounts for 
the variable shape and cross-section of the cometary dust that affects the radiation force 
impulse experienced by the dust. 

We assume the dust size follows a simple power-law with a differential size index of q. 
Therefore, the dust production rate is expressed as 





























iV(rH, ad)dad = No 


—k 




1 /im 


-<? 


dttr. 


(4) 


where A'^o is the mean dust production rate of 1 fim particles and th is the heliocentric 
distance at which th e dust is releas ed. We use k = 4 following the canonical comet 
brightening rate (c.f. 


Everhart 


19671). 


We assume the observed flux is solely contributed by scattered light from the dust 
particles released in the current perihelion passage, and set the start epoch of dust emission 
to 2014 Feb. 18 (r = 50 d) as found in § 12.2.11 Simulated particles are symmetrically 
released around the comet-Sun axis line at the sunlit side. The size distribution is set to 
the interval of the free parameter /5rp,max (i-e. the lower limit of dust size) to an upper size 
limit constrained by the escape speed Uesc = \/‘2GMjq/RG, where Mn = is the 

total mass of the nucleus, u n = 500 kg ■ m~^ th e bulk density of the nucleus, R^ = 1.35 km 


the effective nucleus radius flHowell et ah 


that gas drag be come negligible (iGombosi et ah 


20141). and R g = lOi^N the characteristic distance 


1986). A modihed MERCURY6 package 


(1Chamberslll999l) is used to integrate parti cles from the st art epoch to the observation epoch 


using the 15th order RADAU integrator flEverhart 


19851) . Gravitational perturbations 


from the eight major planets (the Earth-Moon system is represented by a single mass at 
the barycenter of the two bodies), radiation pressure and Poynting-Robertson effect are 
included in the integration. The orbital elements of 209P/LINEAR are extracted from the 
JPL small body database elements 130 (http://ssd.jpl.nasa.gov/sbdb.cgi) and are 
listed in Table |3] together with other parameters used for the model. The resulting image 
is convolved with a 2-dimensional Gaussian function (with FWHM equals to the FWHM 
of the actual images) to mimic observational effects such as the instrumental point spread 
effect and atmospheric seeing. 


We hrst model the May 18 Xingming image with the following procedure. First, we 
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select the pixels > 3a from the background with a the standard deviation of all the pixels 
in the image. Observed and modeled surface brightness prohles are then normalized to 
3 FWHMs beyond the nucleus along the Sun-comet axis, with the region within 1 FWHM 
from the nucleus being masked out, as the signal from the nucleus may contaminate the 
central condensation. The degree of similarity of the two prohles is then evaluated using the 
normalized error variance (NEV), under a polar coordinate system centered at the nucleus 
with angular resolution of 1°: 


, 1 - O, 

NEV = - > ^ ^ 

n ^ 


1=1 


O, 


(5) 


where n is the number of pixels above 3a from the background, Mj and Oj are the pixel 
brightness from the modeled and observed brightness prohle respectively. We set the 
tolerance level of NEV to 5% in order to derive uncertainties of the model parameters. 
We then test a range of parameters as tabulated in Table 01 which yields /9rp,max = 0.005, 
Vo = 40 ± 10 m • s“^, g = 3.8 ± 0.4 and ai^ = 0.3 ± 0.1, shown as Figure |3]). We hnd the 


dominance of larger dust in general agreement with previous results (e.g. 


2014J: 


Younger et ah 


Ye fc Wiegert 


20151), except a steeper size distribution (g = 3.8 vs. g = 3.25) 
and a slightly lower ejection velocity (uej = 1.5 m • s ~^ vs. Vpj = 2 . 5 to 4 .4 m • s“^ for 


millimeter-sized dust) comparing to the results from 


Ishignro et af 


(120151 ). The ejection 


speed is abo ut an order of magnitude lower than the one 


models (e.g. 


Jones 


1995 


Crifo fc Rodionov 


1997 


given 


Williams 


by some classic ejection 


200ll ) and is not much higher 


than the escape velocity (uesc = 0.2 m • s ^). 


We then model the May 25 Gemini image. As the image was taken almost edge-on to 
the comet’s orbital plane, dust particles at different sizes collapse onto the viewing plane, 
making it possible to collapse the image into a 1-dimensional prohle without losing too 
much information. This comes with the beneht of simplifying subsequent analysis. Here we 
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recognize that the orbital plane angle at the time of the observation (+9.3°) was not really 
as small as those nsed in other stndies (generally < 5°), therefore collapsing the image may 
resnlt a loss in the resolntion of the data, which shonld be reflected as an elevation in the 
nncertainties of the modeled parameters. However, we later see that the nncertainties of the 
best models of the May 25 Gemini image are comparable to that of the May 18 Xingming 
image (which was not collapsed). Hence, we conclnde that the collapse of the image does 
not have a signihcant impact to onr resnlt. 

We test the same range of parameters as listed in Table 0] to model the May 25 Gemini 
image. The observed and modeled snrface brightness profiles are then integrated along 
the direction perpendicnlar to the orbital plane, and normalized to 3 FWHMs behind the 
nnclens along the Snn-comet axis. The goodness of the model is determined nsing Eq. |5l 
We hnd /5rp,max = 0.004, Vq = 40 ± 10 m • s“^, q = 3.8 ± 0.4 and = 0.3 ± 0.1, which is in 
good agreement with the parameters fonnd from the May 18 Xingming image. However, 
we note that despite the fact that the £t at the tailward direction is good, the discrepancy 
between the modeled and the observed profile at the snnward direction is striking (Fignre|3]). 
Additional testing at the snnward-only section with the same test grid as Table 0] reveals no 
compatible dnst model (Fignre0]), which snggest a violation of the steady flow assnmption 
we nsed for the model. 


2.2.3. Near Nucleus Environment 


To nnderstand the physical properties of the non-steady coma, we separate the steady 
(i.e. the dnst tail) and the non-steady component (i.e. the coma) in the snrface brightness 


profile and calcnlate the flnx for each of them. We first perf orm an interna. 


photometric calibration nsing the 2MASS stars in the image flSkrntskie et ah 


absolnte 


200d). The 


selected calibration stars are at least 0.5' away from the tail axis to avoid contamination 
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from the comet. We then correlate the observed prohle to the modeled prohle on the 
absolute scale. We subtract the modeled tail component (from 12.2.2p . and interpolate 
the linear portion of the coma to further isolate the nucleus signal (Figure E]). This leaves 
the profiles of the steady tail and the non-steady coma calibrated to an absolute scale. By 
integrating these profiles, we derive the flux for the tail and the coma to be = 0.40 Jy 
and J^coma = 0.04 Jy respectively. The effective cross-section area for each component can 
be calculated by 


C = ( 

A,{a)FQ,^ 


( 6 ) 


where A is the geocentric distance, F\ and Fq^\ is the flux of the component of interest and 
the Sun at the desired wavelength A, which Fq^x = 1.4 x 10^^ Jy for Kg band, and A\{a) 
the phase angle corrected geometric albedo. 


For the tail component, the dust model gives a mean dust size Od = 2 x 10 ^ m. By 


using Aa( 0°) = 0.05 and cal culating the p 


rase angle correction following the compound 


Henyey-Greenstein function fjMarcusll2007al]bl) . we derive Gg = 5 km^ and the corresponding 
dust mass Md = |pdhdGe = 1 x 10® kg. Considering the th dependency, the dust 
production rate at the observation epoch is calculated to be Md = 0.4 kg • s“^, yielding a 

1 : 2 using the water production rate derived fro m narrow banc 


dust-water mass ratio o 

~ 1 

observations ( 

Schleicher 

2OI4I 


2005 


Rotnndi et al. 


This is lower than other measurements (e.g. 


Kiippers et al. 


20151 ) but is perhaps not unexpec ted, given the 


factor of 10-100) of the dust-gas ratio among comets flA’Hearn et al 


arge s catter (within a 


19951) . We also note 


th e derived dust produc tion is about an order of magnitude lower than the value derived 
by 


Ishignro et al. 


(120151 ). likely due to different model parameters (such as a) used for the 


calculation. 


On the other hand, the non-steady region extends no more than ~ 2" behind the 
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nucleus, corresponding to a mean lifetime of ~ 1 d appropriate to 10-100 /im particles. 
Interestingly, this is comparable to the mean l ifetime of icy gra ins (purity Xi > 0.9999) of 


comparable sizes at Th = 1 AU flHanner 


1981 


Beer et al. 


20061. e.g.), seemingly endorsing 


that the presence of an icy grain halo as a self-consistent explanation to the observation. 
However, we note that this hypothesis is not without problems: for dirtier icy grains such 
as Al = 0.9 grains, centimeter-sized grains would be required to survive to 1 d; we also note 
os are known to exist only on long period comets and hyperactive JFCs 


that 

icv grain ha 

os are 

(c.f. 

Gombi et al. 

2013) 


20131) . Therefore, more direct evidence is needed to prove/disprove the 


icy grain halo hypothesis for the case of 209P/LINEAR. 


2.2.4- Nucleus Properties 

As the nucleus is effectively a point source in our data, we reconstruct the nucleus 
signal by htting the isolated nucleus signal in Figure E] with a Gaussian function. This 
yields the nucleus flux Anudeus = 0.02 Jy. As the nucleus size has been reliably measured by 
radar, we derive the corresponding geometric albedo of the nucleus by 


Aa(0°) = 


rn 


-V 


A2 


.lAUy 

where <F(q;) = 10“°'^^“ is the phase angl e function, with a the p hase angle and 
{3 = 0.035 mag • deg“^ the phase slope (e.g. 


(7) 


Gehrels fc Tedesco 


19791) . We yield 


Aa(0°) = 0.12 in the As band. This i mplies a steep spectral slope considering the i?c-band 


albedo constrained by 


As band. This 

mplies 

r Ishiguro et al. 

(20li 


(120151) that is at the order of 0.05, m aking the nuc 


of 20 9P/LINEAR similar to D-type asteroids and most Trojan asteroids flDumas et al 


19981 ). 


eus 
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3. The Meteors 

3.1. Instrument and Data Acquisition 

The Camelopardalid meteor shower was observed using the Canadian Meteor Orbit 
Radar (CMOR). CMOR is an interferometric radar array located near London, Canada. 
The main component of CMOR consists of six stations operated at 29.85 MHz with a pulse 
repetition frequency of 532 Hz. Meteors are detected along a great circle on the sky plane 
perpendicular to the radiant vector, when their ionized trails reflect the radar waves sent by 
the transmitter. Observations are routinely processed by an automatic pipeline to eliminate 


false detecti ons and calculate tr ajectory solu t ions. The c 
be found in 


Jones et al. 

(2005) 

Brown et al. 

to 

o 

o 

00 

) and 


etails of the CMOR op eration can 


Wervk k. Brownl ( 20121) 


In this study, we focus on multi-station data as it allows for reliable determination of 
many meteoroid properties. Single-station data (from the main site) is only used for flux 
calculation. We hrst prepare our initial dataset by extracting Camelopardalid meteors 


from t he processed daily multi-station data, following the procedure described in 


Ye et al 


by 


(120141). The aper t ure (b oth spatial and velocity) are initially set to the predicted value 


Ye fc Wiegert 


(120141) and iterated several times until the optimal va^ 


a ma ximum number of meteors) are found. A Monte Carlo procedure (jWervk &: Brown 


ues (i.e. includes 


20121) is then used to determine the weighted mean radiant and meteor velocity, which are 


found to be A — A© = 38°, /3 = -|-57° in the sun-centered coordinate system and with an 
in-atmosphere velocity Um = 18.8 km • s“^. The sizes of the spatial and velocity apertures 
are then found by comparing to the radiant/velocity density prohle between the outburst 
date and the background as determined from ambient meteor activity ±2 days away from 
the outburst date. As shown in Figure [6l spatial and velocity aperture sizes are determined 
to be 10° and 11% of Um- A total of 99 Camelopardalid meteors are selected in such manner. 
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The meteor population st udied by CMOR can be broadly classified into overdense 


and underdense meteors (e.g. 


McKinley 


19611) . For meteors with similar compositions 


and properties, overdense meteors are typically associated with larger meteoroids and 
vice versa. For the case of the Camelopardalid meteor shower, the size cutoff between 
underdense and overdense meteors is approximately = 0.0003 (equivalent to Od = 2 mm 
assuming pd = 1000 kg ■ m“^). Compared to the underdense meteors, whose appearance are 
usually simple, overdense meteors tend to exhibit a complicated and variable appearance, 
making them sometimes difficult to be identified automatically. Therefore, we retrieved and 
inspected the raw data ±6 hr from the predicted peak of the meteor outburst for overdense 
meteors. A total of 63 Camelopardalid overdense meteors are manually identified in such 
manner, labeled as the overdense dataset. Out of these 63 meteors, 14 of them are also 
found in the initial dataset. We remove these 14 meteors from the initial dataset, leaving 
the other 85 underdense meteors, and label them as the underdense dataset. The three 
datasets are summarized in Table |5l 


3.2. Results and Analysis of the 2014 Outburst 


3.2.1. General Characteristics 


We derive a weighted mean geocentric radiant of oq = 124.9° ± 1.0°, = 79.2° ± 0.2° 

(J2000 epoch) and in-atmosphere velocity Um = 18.8 ± 0.1 km • s“^, using the 99 


Camelopardalid meteors in 


by other studies fj.Tenniskens 


le initial dataset. T 


lis is consistent with the values derived 


2014 : 


Madiedo et ah 


20141 : 


Younger et ah 


20151 ). We also note 


a change in the percentages of overdense and underdense meteors around the peak hour 
(Figure [7]), which may reflect the dynamical delivery of meteoroids at different sizes to the 
Earth’s orbit. 
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We then derive the mass distribution index s (dehned as diV oc m ^dm where m is the 
mass) for the underdense and overdense population respectively. For underdense meteors, 


the cumulative am 

plitude-number relation is tvpicallv used to derive the shower mass index 

(e.g. 

Blaauw et ah 

2011 

); for overdense meteors, 

Te cumulative duration-number relation 

is sometimes used (e.g. 

McIntosh 

1968; 

Ye et ah 

2014 

). For our underdense sample, we 


select 50 underdense meteors with echo range within 110-1 30 km; the range h 


to avoid contamination from overdense transition echoes flBlaauw et ah 


2011 


ter is applied 
I. For the 


overdense sample, all 63 meteors in the overdense d ataset are used. T he data and the 


uncertainty are htted using the MultiNest algorithm flFeroz et al.ll2013l) . taking account the 


number statistics of the data. The technique will be described in a separate paper in more 
detail (Pokorny & Brown, in prep). We hnd Sud = 1-84 ± 0.07 and Sod = 2.02 ± 0.19 for 
underdense and overdense meteors respectively (Figure E]). This can be related to the size 
index q by 


q = 3s — 2 


( 8 ) 


which, for our range of observed s, corresponds to g = 3.5 to 4.1. This agrees with the 
number derived from cometary observations in § 12.2.21 

The flux is calculated from the number of meteors detected per unit time divided 
by the effective collect ing area of the radar system, following the procedure described in 


Brown fc .TonesI fjl995h . The calculation of flux does not require a multi-station setup; 


single-station data is usually sufficient with proper background subtraction. In fact, by 
using the main-station detections, the statistics can be raised by a factor of ~ 5. To derive 
Camelopardalid-only flux, we subtract the raw meteor flux by the background flux following 


the procedure described in 


Ye et ah 


(2013) and ICampbell-Brown fc Brown! (2015). The flux 


is converted to a Zenith Hourly Rate (ZHR) assuming a single power law size distribution 
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applies to the observed size range flKoschack fc Rendtel 


19901). The derived CMOR flux is 
shown in Figure [9] along with the flux derived from visual observation]^. Overall, radar and 
visual observations show agreement in terms of activity timing, with a moderate rise and a 
steep decline in rates, as well as a main peak around 8h UT, 2014 May 24. We note that 
the visual profile suffers from small statistics (only ~ 15 meteors per bin during the peak, 
comparing to ~ 60 meteors for the radar), and so the two “peak-lets” at 6:30 and 8:30 UT 
are likely artifacts. In both techniques, further refinement of the exact peak time is perhaps 
not meaningful due to the relatively small statistics of the data. The CMOR flux (corrected 
to a limiting magnitude of -1-6.5) is about half an order of magnitude higher than the visual 
flux, seemingly indicating an overabundance of faint meteors and a break in the power law 
somewhere beyond the naked-eye limit. 


3.2.2. Meteoroid Properties 


The Camelopardalids have almost identical entry speeds and geometry (with respect 
to CMOR) as another JFC shower, t he October Draconid s, which was observed by CMOR 


during its 2011 and 2012 outbursts (lYe et ah 


2013 . 


2014J) . This coincidence allows us to 


directly compare the main characteristics of these two showers independent of instrumental 
effects or entry speed corrections. A distinct difference between the two showers is in the 
specular height distribution of the meteors: the Draconids appear 5-10 km higher than 
the Camelopardalids as observed by CMOR (Figure fTOl). It has long been thought that 


the exceptional ablation heie 

it of 

;he Draconids is the direct consequence of their extreme 

). Hence, a simple interpretation of the observed height 

fraeilitv fe.e. Borovicka et ah 

2007 



^Available at http://www.imo.net/live/cameleopardalids2014/, retrieved on 2015 


April 2. 
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distribution of the two showers is that the Camelopardalid meteoroids are less fragile 
relative to the Draconids. As the outbursts from the two showers originated from cometary 
ejecta with young ejection ages (less than a few hundred years), the difference in space 
weathering is not significant; the observations seem to suggest that the surface material 
properties of the two parent comets are different. 


We compare our result to the results derived from other Camelopardalid studies. 


Younger et al. 


(120151) . who also observed the 2014 Camelopardalid outburst with a meteor 


radar, reported that the Camelopardalid meteoroids were less fragile than sporadic 
meteoroids, a hnding that is not apparent in our Figure [TUI due to our aggressive binning 
to enhance the statistics; but Younger et ah’s finding is at least qualitatively consistent 


with our finding that the Camelopardalid me teoroids being less fra gile relative to 


Draconids. Conversely, optical observations by 


■Tenni.skensl (120141 ) and 


Madiedo et al. 


he 


(12014 


show that the Camelopardalid meteoroids are very fragile and are consistent with fluffy 
aggregates like the Draconids. However, we note that (1) optical observations are sampling 
meteoroids of a larger size range (close to centi meter-sized , while rada r observations a re 


samp 


ing millimeter-sized meteoroids); and (2) Ijenniskensi (12014 and 


Madiedo et al. 


(12014 ’s observed meteors were recorded in a wider time span than the radar (on the order 


of 1 d vs. a few hours). Meteors detected away from the predicted peak mainly consist of 
background meteoroids that are part of older, disrupted trails. Hence, the optical meteors, 
whose properties seem very different than the radar meteors, may represent Camelopardalid 
meteoroids at different sizes and ages. 


3.3. Camelopardalid Activity in Other Years 

We conduct a search in the CMOR database for any undetected Camelopardalid 
activity in previous years, using the 3-dimensional wavelet analysis technique (e.g. 
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Brown et al. 


2010 


Bruzzone et al. 


2 OI 5 I) to compute the wavelet coefficient at the location 


of the Camelopardalid radiant. The time window is restricted to one week around the nodal 
passage of 209P, namely in the solar longitude range Aq = 60° — 66°. CMOR has been fully 
operational since 2002, but data in 2002, 2006, 2009 and 2010 are severely (off-line periods 
more than 24 hours) interrupted by instrumental issues; hence we only inspect years with 
complete data for possible Camelopardalid activity. 


We hnd distinct activity in 2011, while the activity in other years, if any, was too 
weak to be reliably separated from the background (Figure [TT]). The 2011 outburst is even 
noticeable on the raw, unprocessed radiant map fFigure [T^ . albeit much weaker than the 
2014 outburst. We were able to extract 15 meteors for the 2011 outburst, which yields a 
weighted radiant of = 119.5° ± 2.1°, Aq = 77.2° ± 0.3° (J2000 epoch) and in-atmosphere 
velocity Um = 19.3 ± 0.3 km • s“^. We hnd no obvious peak of activity, but the core of the 
activity falls between 2011 May 25 at 6-11 h UT (Aq = 63.6°). The 2011 activity was not 
high enough to derive a statistically meaningful hux, but we estimate the 2011 hux to be 
about an order of magnitude lower than the 2014 hux, since the number of raw echoes is 
roughly 1/10 of that of 20li^. By following the same technique described in § 13.2.11 we 
derive a la upper limit of the hux to be < 0.01 km“^ • hr“^ for other years. 


3.4. Modeling the Dust (II) 

The dust model derived from cometary observations has placed some useful constraints 
on the physical properties of the Camelopardalid meteoroids. In this section, we explore 
the contribution of young meteoroid trails (dehned as trails formed within ~ 50 orbital 

^The change of radar collecting area in diherent years is negligible thanks to the high 
declination of the Camelopardalid radiant. 
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revolutions) to the observed meteor activity using numerical techniques. Older dust trails 
have experienced more perturbations from the major planets and are too disrupted to 
model. The simulation procedure is essentially the same as that in § 12 .2.21 apart from 
extending the integration time several hundred years backward. To address possible meteor 
act ivity, we select a subse t of Earth-approaching mete oroids following the method discussed 


by 


Brown fc JonesI (1199811 and 


Vaubaillon et ah 


(I2nn5h : 


AX = Urel X AT 


(9) 


where Urei ~ 17 km • s“^ is the relative velocity between the meteoroid and the Earth, 
AT is the characteristic duration of the meteor shower which we take as AT = 1 d. These 
yield AX = 0.01 AU. The simulated meteoroid is included in the subset when its Minimum 
Orbit Interse ction Dis t ance fMOID) to the Earth’s orbit, calculated with the subroutine 


developed by 


Gronchil (120051) . is smaller than AX. 


We use the dust model derived from o ur cometary observ a tions for the ejection of 


meteoroids. For comparison, the traditional 


Crifo fc RodionovI (119971) model (denoted as 


the C&R model hereafter) is also used in a parallel simulation. The start of the integration 
is set to 50 orbits ago (or about 1750 A.D.). We hrst integrate 209P/LINEAR back to 
the year of 1750, and then integrate it forward with meteoroids released at each perihelion 
passage when the parent has rn < 1.4 AU, the heliocentric limit of cometary activity as 
indicated by cometary observations. When the simulation is hnished, we examine the 
encounters of all meteoroid trails in the years that CMOR was operational. 


The results from both ejection models are largely identical, making it difficult to 
distinguish the better ejection model using observations. This emphasizes that the 
evolution of older trails is predominantly controlled by planetary perturbations rather than 
ejection speed. The 2014 encounter is easily identihable thanks to the high density of the 
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corresponding trail (Figure [13]), with the simulation agreeing with the observations. We 
also note that our simulation predicts the Earth would hrst encounter larger meteoroids 
1 Figure [TB]1 . a result consistent with CMOR observation of early overdense meteors noted 
in ^ l3Xn 

The flux of meteoroids can be estimated by relating the number of meteoroids in 
Earth’s vicinity to the dust production rate of the comet. From the analysis in § 12.2.31 
we estimate the current dust production rate of 209P/LINEAR is of the order of 10® kg, 
or ~ 10^^ meteoroids per orbit (taking a ~ 10“^ m as found previously). From the 
meteoroid stream simulation, we hnd ~ 1% of the meteoroids released between 1750-2014 
are delivered to the Earth’s vicinity during the 2014 encounter, corresponding to a flux of 
T^l%x N X AX-2 ^ ^ ^ ^q 14 ^ X (1 d)-^ = 0.02 km-2 ■ hr’^ 

comparable to the flux determined from visual and radar meteor observations. This implies 
that 209P/LINEAR was not substantially much more active in the past several centuries, 
an idea also supported by the apparent lack of annual activity of the Camelopardalid 
meteor shower. 

Additionally, we hnd predicted encounters in 2004, 2008 and 2011 from our simulations 
(Figure mi). The 2004 and 2008 encounters are predicted to be about an order of magnitude 
weaker than the 2011 encounter, thus we expect this activity to be buried in the sporadic 
background. The 2011 case is interesting as the parent was near aphelion at the time 
of the meteor outburst. Both the C&R model and our ejection model derived from the 
cometary observations only indicate encounters with a few extremely weak trails formed 
between 1763-1768 in 2011. The calculated peak time and width (both ejection models 
suggest peak times of 2011 May 25 ~5:40 and 9:00 UT for the 1763- and 1768-trail, with 
full-width-half-maximum of ~ 8 hr) is consistent with CMOR observations. However, the 
hux predicted by the model is by a factor of 100 lower than what was observed, potentially 
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hinting at a significant but transient increase of activity of 209P/LINEAR around those 
epochs. The same 1763- and 1768-trail also contribute to the 2014 meteor event; however, 
the overlapping peak time between trails (mostly < 1 hr apart) makes it difficult to 
distinguish activity from individual trails in the observations. 


Discussion 


4.1. The Dynamical Evolution of 209P/LINEAR 


Recent work by iFernandez &: Sosal (1201411 revealed a set of unique members among the 
JFCs that reside in highly stable (> lO'^ yr) orbits, including 209P/LINFAR. We extend 
their work for the case of 209P/LINFAR by generating 1000 clones of 209P/LINFAR 
using the orbital covariance matrix provided in JPL 130, and integrate all of them 10® yr 


backwards. The in 


integrator fiRnlirsch 


; egrat i on is perfor med with MERCURY6 using the Bulirsch-Stoer 


1972 


Stoer 


197211 . 


As shown in Figure [TSl the core of the clones remain in Earth’s vicinity for ~ 10^ years, 
much longe r than the typical p hysical lifetime for similar-sized JFCs in the near-Earth 


region (e.g. 


DiSistoetah 


200911 . In addition, we note the core of the clones is extremely 


compact for more than 100 orbits (la width in semimajor axis ~ 0.0002 AU), until an 
extreme close approach to the Earth (dmin ~ 0.006lo;oo5 AU) around 1400 Mar 12 (on 
Julian calendar) scatters the clones. The miss distance of this approach to t he Earth is 


smaller than the recorded close approach by Lexell’s Comet in 1770 (0.015 AU; 


Kron 


200811 


and prompted us to look at medieval astronomical records for possible sightings, without 
success. If the activity level of 209P/LINEAR in the 15th century is comparable to what it 
is now, the comet would have been -|-7 mag during its approach in 1400, below the naked-eye 
limit of medieval astronomers; however, any significant (by several magnitudes) increase in 
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activity could have been noticeable. The lack of possible sightings for 209P/LINEAR’s close 
approach in 1400 suggests that the comet was not substantially more active ~ 100 orbits 
ago. 


Since 209P/LINEAR is in a stable orbit, the associated meteoroid stream may also 
possess a set of orbits that are more stable than other JFC streams. To quantify the 
dispersion process of the Camelopardalid meteoroid stream, we adopt the same integration 
procedure as described in § 13.41 and examine the evolution of meteoroid trails released 
between 1-revolution (5 yr) and 1000-revolution (5000 yr), shown as Figure ITTl It can be 
seen that the narrow stream structure is maintained for trails that formed as far as ~ 1000 
to 2000 yr ago, which i s a few times longer th an other JFC streams such as the vr-Puppid 


meteoroid stream (e.g. 


Cremonese et ah 


19971) . We also note that the meteoroid stream 


evolves differently than the parent. The degree of the difference increases as the age of the 
stream increase. For example, the current radiant of the core of 200-rev meteoroids (i.e. 
meteoroids released at about 1000 A.D.) would be at oq = 120°, Jq = +60°, encountered 
at Aq = 70° (approximately June 1). There is no established meteor activity related to this 


hypothetical radiant, alth ough a few other possib 


with 209P/LINEAR (e.g. 


Rudawska &: Jenniskens 


e annual showers 

rave 

2014; 

Segon et ah 

2014) 


Deen associated 


4.2. Nature of 209P/LINEAR and Comparison with Other Low Activity 

Comets 

Following our analysis, it seems evident that 209P/LINEAR has been mostly weakly 
active for the last few hundred orbits, while it might have been in a near-Earth JFC orbit 
on the time scale of 10^ yr. This is compatible with the idea of 209P/LINEAR as an 
aging comet exhausting its remaining near surface volatiles as derived from the classical 
interpretation of cometary evolution. It is perhaps not possible to know how long the comet 
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has stayed in the inner solar system; however, we note that the gradual decrease of the 
perihelion over the course of few thousand years (as indicated in Figure USD may provide a 
prolonged favorable environment for weak cometary activity, as the sub-surface volatiles 
undernea th the dust mantles c an be (re-)activated by the gentle decrease of the perihelion 


distance flR.ickman et al. 


199nh 


What does 209P/LINEAR tell us about other low activity comets? In the following we 
briefly discuss three other Earth-approaching comets (i.e. those that may generate meteor 
showers) listed in Table [T] and compare them to 209P/LINEAR. The other hve comets in 
the list do not generate meteor showers, making it difficult to address their physical history 
in a manner similar to 209P/LINEAR. 


252P/LINEAR Little is known about this newly discovered comet at the moment, ex- 
cept that nume rical simulation indicate a recent (< 100 orbits) entry to the inner solar system 


fjTancredill20l4 see also http: //www. astronomia. edu.uy/Cr iter ion/Comet s/Dynamics/table_num.htm: 


retrieved 2015 May 17), implying a different origin and evolution compared to 
209P/LINEAR. Considering its young dynamical age in the inner solar system, the low 
activity of 252P/LINEAR may reflect a relative lack of volatiles at the time of formation of 
the nucleus. 


289P/Blanpain 289P/Blanpain is the only low acti vity comet in th e list that is 


associated unambiguously with annual meteor activity (jjenniskens 


2008). The comet itself 


was lost for some 200 yr after its initial discovery in 1819 (and had been referr ed as D/1819 


Wl), until being re-discovered as a faint asteroidal body 2003 WY25 in 2005 flFoglia et al. 


20051 ). Multiple clues suggest 2003 WY25 is the r emnant of the original 289P/Blanpai n 


follow ing a catastrophic fragmentation event (e.g. 


■lenni.skens fc Lyytinen 


2005 


Jewitt 


20061 ). Hence, the low activity nature of 289P/Blanpain may have a completely different 
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origin than that of 209P/LINEAR. 


300P/Catalina 300P/Catalina (known as 2005 JQ5 in some early literatures) is 
interesting, as it is the only othe r come t in our list that is concurrently classihed as a 


stable JFC by 


Fernandez fc Sosal (120141) . It has not been associated with any established 


annual meteor shower, althou gh a few possible linkag es have been sugges 


Rudawska fc Jenniskens 


2014() . Radar observations by 


Harmon et ah 


ed (e.g. 


( 20061 1 revealed a 


rough surface similar to 209P/LINEAR; however, the presence of cm-sized c 


nucleus of 300P/Catalina, which is absent for 209P/LINEAR ( Howell et al 


ust around the 


201411 . seems 


to indicate stronger outgassing activity of 300P/Catalina compared to 209P/LINEAR at 
the present time. It may be possible that 300P/Catalina is at an earlier stage of dormancy 
compared to 209P/LINEAR. 


5. Conclusions and Summary 

The low activity comet, 209P/LINEAR, may indeed be an aging comet that is quietly 
exhausting its last bit of near surface volatiles. This idea is supported by the convergence 
of several different lines of evidence: dust modeling of cometary images that revealed a 
presently weakly active comet, analysis and modeling of meteor observations that revealed a 
low dust production over the past few hundred orbits, numerical analysis of the dynamical 
evolution of the comet that suggested a stable orbit in the inner solar system over a time 
scale of 10'^ yr. 

The main Endings of this paper are: 

1. The best-fit dust model to the cometary images involves a low ejection speed (1/10 of 
moderately active comets) and large dust grains (fid = 10“*^ m). The dust production 
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rate of the comet at 19 d after perihelion is 0.4 kg • s“^, a remarkably small number. 

2. The coma region appears to be inconsistent with the steady-flow model. The general 
characteristics of this region is compatible with the icy grain halo theory, a theory 
that is known to be only applicable to active long period comets and hyperactive 
Jupiter-family comets. More conclusive evidence is needed to establish or disprove 
this hypothesis. 

3. By applying a coma subtraction technique, the nucleus signal is separated from the 
coma, yielding a geometric albedo Aa( 0°) = 0.12 appropriated to band. Coupling 
with optical measurements at visible band, this indicates a reddish spectrum of the 
nucleus of 209P/LINEAR similar to that of D-type asteroids and most Trojans. 

4. Radar observations by CMOR show the peak of 2014 Camelopardalid meteor outburst 
around 2014 May 24 at 8 h UT. From CMOR observations, we derive a mean radiant 
of ttG = 124.9° ± 1.0°, Jg = 79.2° ± 0.2° (J2000 epoch), mean in-atmosphere velocity 
Urn = 18.8 ± 0.1 km ■ s“^, and a peak flux of 0.06 km“^ • hr“^, consistent with visual, 
optical and other radar observations. Numerical simulation conhrms that the outburst 
originated from the dust trails formed in the 18-20th century, a time that the parent 
was perhaps not much more active. The mass distribution index of the meteors, 

s = 1.8 to 2.0, agrees the size index q = 3.8 derived from the modeling of the cometary 
images. 

5. A direct comparison to the Draconids, a meteor shower with almost identical entry 
speed that was also observed with CMOR, shows that a distinctly different height 
distribution between the Camelopardalids and Draconids: the Camelopardalids tend 
to appear ~ 10% lower than the Draconids. This is likely due to the Camelopardalids 
being less fragile relative to the Draconids, the latter of which have long been known 
as extremely fragile meteoroids. This agrees with other radar measurements but 
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differs from optical measurements, which support highly fragile meteoroids. As optical 
observations are sampling meteoroids at larger sizes and wider arrival times, the 
difference in meteoroid properties derived from different techniques may be due to 
sampled meteoroids of different sizes and ages. 

6. We examine CMOR data from 2003 onwards (except 2006, 2009 and 2010) and hnd 
a previously unnoticed Camelopardalid outburst in 2011. The activity peaks around 
2011 May 25 between 6-11 h UT, with a peak flux of the order of 0.005 km“^ ■ hr“^. 
Numerical simulations suggest the dust trail encountered in 2011 was formed in 
1763-1768, however the predicted flux seems to be by a factor of 100 smaller than 
what was observed. This may indicate some temporary increase in activity of 
209P/LINEAR around those times. 

7. Numerical integration indicates 209P/LINEAR may have resided in a stable 
near-Earth JFC orbit for ~ 10^ yr. The dispersion time scale for the Camelopardalid 
stream is about 1000-2000 yr, which is a few times longer than JFC streams such as 
the TT-Puppids. The lack of signihcant annual activity of the Camelopardalid shower 
may serve as a strong evidence of the low activity of 209P/LINEAR over the past 
several hundred orbits. 

8. We compare 209P/LINEAR to three other low activity comets that are associated 
with known or hypothetical meteor showers: 252P/LINEAR (associated with a 
hypothetical meteor shower in the constellation of Lepus), 289P/Blanpain (associated 
with the Phoenicid meteor shower), and 300P/Catalina (associated with a few 
possible meteor showers, such as the June e-Ophiuchids). A diversity is seen: the 
low activity of 252P/LINEAR may be congenital; that of 289P/Blanpain may be 
due to catastrophic fragmentation. 300P/Catalina shares many similar physical 
and dynamical characteristics with 209P/LINEAR; but the presence of cm-sized 



meteoroids around the nucleus may indicate a stronger outgassing activity of 
300P/Catalina compared to 209P/LINEAR at the moment. 
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+ 4 " + 2 " 0 " 


R.A. to Nucleus 

Fig. 1.— The 2014 Apr. 9 GMOS-N image (stretched in logarithm scale) superimposed with 
the synchrone model. The ages of the synchrone lines (dashed lines) are (in counterclockwise 
order) 10, 25, 50 and 100 d respectively. The oldest visible dust was released at r ~ 50 d, 
appropriate to late Feb. 2014. 
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Fig. 2.— Composite images of 209P/LINEAR taken by Xingming 0.35-m telescope and 
Gemini Flamingo-2 on 2014 May 18 and 25. The images are stretched in asinh scale and are 
rotated to have north-up east-left. 
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1 " 10 " 

Distance to Nucleus 


Fig. 3.— Observed (colored pixels) and modeled (contours) surface brightness profiles for the 
Xingming image (upper hgure) and the Gemini F-2 image (lower hgure; the sunward data 
is shifted downwards for clarity). The surface brightness prohles are normalized to the pixel 
intensity 3 FWHMs behind the nucleus along the Sun-comet axis to avoid contamination 
from the nucleus signal. The mean best model for both the Xingming and the Gemini F-2 
images has /drp.max = 0.004 to 0.005, Vq = 40 m • s“^, q = 3.8 and = 0.3. 
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Fig. 4.— Representative attempts to fit the sunward section of the coma in the Gemini F-2 
image. The observed and modeled profiles are all normalized to 3 FWHMs away from the 
nucleus along the comet-Sun axis. These models have q = —3.8 and /5rp,max = 0.004. 
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Fig. 5.— Separation of the coma and nucleus signal based on the Gemini F-2 image. Upper 
hgure: observed prohle and modeled prohle from the dust model. Middle hgure: derived 
coma+nucleus prohle by subtracting the observed prohle with the modeled prohle. Lower 
hgure: nucleus-only prohle, derived from subtracting the linear portion of the coma prohle. 
The X-axis corresponds to the Sun-comet axis. 
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Velocity Aperture Size (%) 

Fig. 6.— Determination of the optimal radiant and velocity apertures. Radiant aperture 
is centered at A — Aq = 38°, (3 = +57° in the Sun-centered ecliptic coordinate system, 
(in-atmosphere) velocity aperture is centered at Um = 18.8 km ■ s“^. Background values are 
extracted from non-outburst dates ±2 days from the outburst date (i.e. 2014 May 22 and 
26). The optimal radiant and speed apertures are determined to be 10° and 11% respectively 
(marked by arrows). The velocity aperture is determined for the spatial aperture of 10°. 
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Fig. 7.— Top: Variations of the overdense meteor fraction with Poisson errors, binned in 
2 h intervals. A dip (i.e. larger proportion of small meteoroids) is apparent aronnd the 
peak honr (7-8h UT). Bottom: Raw nnmbers of overdense and nnderdense Camelopardalid 
meteors detected by CMOR, binned in 15 min intervals. 
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Fig. 8.— Determination of mass indices for the underdense (upper figure) and overdense 
(lower figure) populations. The mass indices are determined to be 1.84 ±0.07 for unserdense 
and 2.02±0.19 for overdense meteors. The dashed lines show the best fit as determined by the 
technique developed by (Pokorny & Brown 2015, in prep). The uncertainties are based on the 


distri butions of the posterior probabilities obtained by the Mult i Nest a lgorithm flFeroz et al. 
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Fig. 9.— The variation of the flux (corrected to a limiting magnitude of +6.5) of the 2014 
Camelopardalid meteor outburst as observed by CMOR and IMO visual observers. The 
CMOR observations are binned in 1 hr intervals. Error bars denoting Poisson errors. 
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Fig. 10.— Specular height distribution of the underdense meteor echoes observed by 
CMOR for the 2011/12 Draconid outbursts (denoted as DRAll and DRA12) and 2014 
Camelopardalid outburst (denoted as CAM14), plotted as shaded bars. Specular height 
distribution of sporadic meteors (generated using all meteors detected by CMOR with 
within 5% from 20 km • s“^) is shown as line. 
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Fig. 11.— Variation of the relative wavelet coefficient at A — Aq = 38°, /3 = +57° and 
u = 20 km • within Aq = 30° — 90° in 2003-2014 (except 2006, 2009 and 2010). The 
expected Camelopardalid activity period is shaded. Activity is noticeable only in 2011 and 
2014. 
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Solar Longitude 


Fig. 12.— Upper figure: the raw radiant map of all meteor echoes detected by CMOR on 
2011 May 25, corresponding to solar longitude Aq = 63°. Angular axis represents R.A. and 
the radial axis represents Declination, both in geocentric coordinates in J2000 coordinates. 
Radiants are plotted as black dots. The Camelopardalid activity is clearly visible near ac = 
120°, Ag = +80°. Lower figure: variation of the relative wavelet coefficient at A — A© = 38°, 
(3 = +57° and w = 20 km • s“^ in 2011, with the Camelopardalid activity marked by an 
arrow. Solid and dashed lines are median and 3(T above median, respectively. 
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Fig. 13.— Nodal footprint of the 1750-2000 trails around 2014 May 24, using the e j ection 


model derived from comet observations (upper hgure) and the ICrifo fc R.odionovl (119971 ) 
ejection model (lower hgure). 
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Fig. 14.— Nodal footprint of the 1750-2000 trails around 2004 May 24, 2008 May 25 and 


201 

the 


May 25, using t h e ejec tion model derived from comet observations (upper row) and 


Crifo fc Rodionov! (119971 ) ejection model (lower row). The scale of meteoroid number is 


identical to that of Figure IT^ but for clarity the meteoroids in this hgure are marked with 
larger symbols. 
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Fig. 15.— Dynamical evolution of 1000 clones of 209P/LINEAR in a time interval of 10^ yr 
with a zoomed section for within 1000 yr. The median (black line) and ila region (shaded 
area) is shown. A highly stable section is seen up to 3 x 10^ years, of which the core of the 
clones remain in near-Earth region and 95% of the clones remain in bounded orbits. 
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Fig. 16.— The arrival distribution of large, overdense-like (od = 5 mm or (],.p = 0.0001) and 
small, underdense-like fan = 1 mm or drn = 0 .0005) meteoroids from observation-derived 


(upper figure) and the 


Crifo &: Rodionovi fjl997l) ejection models (lower hgure) for the 2014 


Camelopardalid meteor outburst. It is apparent that larger meteoroids arrived earlier than 
smaller meteoroids, consistent with CMOR observations. 
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Fig. 17.— Evolution of secular orbital elements of meteoroids of different ages: 1-rev (mete¬ 
oroids released 5 yr ago), 40-rev (released 200 yr ago), 100-rev (released 500 yr ago), 200-rev 
(released 1000 yr ago), 400-rev (released 2000 yr ago) and 1000-rev (released 5000 yr ago). 
The meteoroid ejection model is based on comet observations, but the result is insensitive 
to the choice of ejection model, as the evolution of meteoroid stream is predominantly con¬ 
trolled by planetary perturbations over the investigated time scale. It can be seen that the 
dispersion time scale of the Camelopardalid meteoroid stream is at the order of 1000 yr 
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Table 1: A list of low activity comets according to the definition given in §1. 


Comet 

Ml 

i?N 

(km) 

Assoc, meteor shower 

lOP/Tempel 2 

13.2 

10.6'^ 

- 

28P/Neujmin 1 

11.5 

21.4'" 

- 

102P/Shoemaker 1 

15.7 

1.6^" 

- 

184P/Lovas 2 

14.4 

6.2^ 

- 

209P/LINEAR 

16.9 

2.7^= 

Camelopardalids 

252P/LINEAR 

18.6 

0.5'^ 

Predicted, not yet observed® 

289P/Blanpain 

22.9 

0.32^^ 

Phoenicids 

300P/Catalina 

18.3 

1.4^ 

June e-Ophiuchids (?) 

C/2001 OG108 (LONEOS) 

13.1 

13.6'" 

- 


“ Lamv et al. 

2004 

. 

^ Scotti 

1—1 

. 

'Howell et al. 

(20U 

!). 


'^Drahus (2015, personal communication). 


' Jewitt 

to 

o 

o 

^ _ 

)■ 

t Harmon et al. 

(2006) 


^Unpublished data from Maslov (http : //f eraj .narod.ru/Radiants/Predictions/252p-ids2016eng.html, 
retrieved 2015 May 2). 
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Table 2: Summary of the imaging observations of 209P/LINEAR. 


Time (UT) 

Facility 

Res. 

km/pix 

Exposure 

min. 

Airmass 

FWHM 

arcsec 

ru 

AU 

Z\ 

AU 

Plane Angle 

2014 Apr 9.25 

GMOS-N 

23 

0.2 

1.7 

0.4 

1.043 

0.441 

-34.7° 

2014 May 18.75 

XM 0.35-m 

77 

84 

1.3 

4.4 

0.986 

0.117 

-16.9° 

2014 May 25.94 

F-2 

8 

10 

1.7 

0.8 

1.009 

0.064 

+9.3° 
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Table 3: Input parameters for the Monte Carlo dust model. The orbital elements are ex¬ 
tracted from the JPL elements 130, epoch 2011 Jun 8.0 UT. 


Parameter 

Value 

Semimajor axis a 

2.93102 AU 

Eccentricity e 

0.69237 

Inclination i 

19.44783° 

Longitude of the ascending node 

65.46431° 

Argument of perihelion uj 

150.46931° 

Epoch of perihelion passage tp 

2009 Apr 17.43973 UT 

Nucleus radius i?N 

1.35 km 

Nucleus bulk density pn 

500 kg ■ m“^ 

Dust bulk density pd 

1000 kg ■ m“^ 

Dust albedo 

0.05 
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Table 4: Dust model parameters derived from observations of Xingming 0.35-m (XM) and 
Gemini F-2 (F-2). 


Parameter 

Tested Values 

Best-fit Values 

Dust size lower limit, /9rp,max 

0.001-0.1 in steps of 

XM: 0.005 


1/40 of full range in log space 

F-2: 0.004 

Mean speed of = 1 dust at 1 AU, Vq 

10-400 in steps of 10 

XM & F-2: 40 ± 10 m ■ s"! 

Lagging parameter, 

0.0-0.5 in steps of 0.1 

XM & F-2: 0.3 ±0.1 

Size index, q 

2.6 to 4.4 with steps of 0.1 

XM & F-2: 3.8 ±0.4 
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Table 5: Summary of the CMOR datasets used for analyzing the 2014 Camelopardalid out¬ 
burst,_ 

Label N Description 

Initial 99 Extracted from processed daily data, 

includes 85 underdense meteors and 14 overdense meteors. 
Underdense 85 Subset of initial dataset, contains only underdense meteors. 
Overdense 63 Manually extracted from raw data, 

includes 14 meteors from initial dataset. 



